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Abstract 


In  this  memo,  the  probability  density  functions  (pdfs)  for  the  phase  and  the  envelope  of 
an  amplitude-modulated  (AM)  signal,  where  the  modulation  is  modeled  as  a  Gaussian 
random  process,  in  the  presence  of  narrow-band  Gaussian  noise  arc  calculated.  An 
exact,  closed-form  expression  is  obtained  for  the  phase  pdf,  while  approximate  results, 
accurate  for  small  and  large  carrier-to-noise  ratios  (CNRs),  arc  derived  for  the  envelope 
pdf. 


Resume 


Dans  le  present  document,  on  etablit  les  fonctions  de  densite  de  probabilite  (fdp)  de  la 
phase  et  de  l’enveloppe  d’un  signal  module  en  amplitude  (AM),  dont  la  modulation  est 
modelisee  sous  forme  d’une  fonction  aleatoire  gaussienne,  en  presence  de  bruit 
gaussien  a  bande  etroite.  Une  expression  analytique  exacte  est  obtenue  pour  la  fdp  de  la 
phase,  alors  que  des  resultats  approximatifs,  precis  pour  des  rapports  porteuse/bruit 
faibles  et  eleves,  sont  obtenus  pour  la  fdp  de  l’enveloppe. 


DRDC  Ottawa  TM  2004-037 


I 


This  page  intentionally  left  blank. 


DRDC  Ottawa  TM  2004-037 


Executive  summary 


Amplitude  modulation  (AM)  is  one  of  the  most  common  modulation  types  in  civilian 
and  military  communications  applications.  It  plays  an  important  role  in  electronic 
support  measures  (ESM),  electronic  intelligence  (ELINT)  and  communication 
intelligence  (COMINT).  In  this  memo,  the  probability  density  functions  (pdfs)  of  the 
envelope  and  the  phase  of  AM  signals,  where  the  modulation  is  assumed  to  be  a 
Gaussian  process,  arc  investigated.  The  derived  results  can  be  used  in  the  following 
applications.  The  first  is  the  identification  of  intercepted  AM  communication  signals 
based  on  the  statistical  properties  of  their  envelope  and  phase.  The  second  is  the 
modeling  of  intercepted  radar  or  sonar  signals  as  having  random  amplitudes  even 
though  the  transmitted  signals  arc  of  constant  amplitude.  This  random  amplitude 
accounts  for  “multiplicative  noise,”  caused  by,  for  example,  dispersion  in  the 
propagation  medium.  In  these  two  applications,  the  statistical  properties  of  the 
envelope  and  the  phase  of  the  signals  have  not  been  well  treated  due  to  a  lack  of  the 
appropriate  pdfs.  In  this  study,  a  closed-form  pdf  for  the  phase  is  derived  as  a  function 
of  the  carrier-to-noise  ratio  (CNR)  and  the  signal  parameters.  For  the  envelope  pdf, 
approximate  expressions  arc  obtained  for  small  and  large  CNRs.  These  approximate 
pdfs  arc  compared  with  the  numerically-integrated  results. 
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Sommaire 


La  modulation  d’ amplitude  (AM)  est  un  des  types  de  modulation  les  plus  couranmient 
utilises  dans  les  applications  des  communications  civiles  et  militaires.  Elle  joue  un  rle 
important  dans  les  mesures  de  soutien  electronique  (MSE),  le  renseignement 
electronique  (ELINT)  et  le  renseignement  sur  les  communications  (COMINT).  Dans  le 
present  document,  on  etudie  les  fonctions  de  densite  de  probabilite  (fdp)  de  l’enveloppe 
et  de  la  phase  de  signaux  AM  pour  lesquels  on  suppose  que  la  modulation  est  un 
processus  gaussien.  Les  resultats  obtenus  peuvent  etre  utilises  dans  les  applications 
mentionnees  ci-apres.  La  premiere  consiste  a  identifier,  en  se  basant  sur  les  proprietes 
statistiques  de  leur  enveloppe  et  de  leur  phase,  des  signaux  de  communication  AM 
interceptes.  La  deuxieme  consiste  a  modeliser,  sous  forme  de  signaux  d’ amplitude 
aleatoire,  des  signaux  radar  ou  sonar  interceptes  qui  ont  neanmoins  ete  emis  avec  une 
amplitude  constante.  Cette  amplitude  aleatoire  tient  compte  du  “bruit  multiplicatif  ’ 
cause,  par  exemple,  par  la  dispersion  dans  le  milieu  de  propagation.  Dans  ces  deux 
applications,  les  proprietes  statistiques  de  l’enveloppe  et  de  la  phase  des  signaux  n’ont 
pas  ete  traitees  correctement,  en  raison  de  T  absence  de  fdp  appropriees.  Dans  la 
presente  etude,  on  etablit  une  fdp  analytique  pour  la  phase  en  fonction  du  rapport 
porteuse/bruit  et  des  parametres  des  signaux.  Dans  le  cas  de  la  fdp  de  l’enveloppe,  des 
expressions  approximatives  sont  obtenues  pour  des  rapports  porteuse/bruit  faibles  et 
eleves.  Ces  fdp  approximatives  sont  comparees  aux  resultats  obtenus  par  integration 
numerique. 
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1.  INTRODUCTION 


Amplitude  modulation  is  one  of  the  most  common  modulation  types  in  civilian  and 
military  communications  applications.  It  also  plays  an  important  role  in  electronic 
support  measures  (ESM),  electronic  intelligence  (ELINT)  and  communication 
intelligence  (COMINT)  [2].  In  this  memo,  the  pdfs  of  the  envelope  and  phase  of  AM 
signals  are  investigated,  where  the  modulation  is  assumed  to  be  Gaussian.  The  results 
can  be  useful  in  the  following  two  applications.  The  first  application  is  the  exploration 
of  statistical  properties  of  envelope  and  phase  of  AM  signals  for  signal  classification 
and  identification  [3],  The  second  application  is  for  modeling  radar  and  sonar  signals  as 
having  random  amplitudes  even  though  the  original  signal  is  constant  amplitude  [4]  [5]. 
The  random  amplitude  is  caused  by  dispersion  in  the  medium  and  is  often  referred  to  as 
“multiplcative  noise”.  The  random  amplitude  model  is  considered  to  be  an  appropriate 
model  in  radar-  and  sonar  applications  because  it  represents  the  real  world  [6],  However, 
because  of  the  lack  of  appropriate  representation  for  the  pdfs  of  the  envelope  and  phase 
of  an  AM  in  the  literature,  statistical  properties  of  the  envelope  and  phase  of  an  AM 
have  not  been  well  treated  or  just  ignored.  An  exact  solution  appeal's  to  be  difficult. 
However,  a  good  approximation  would  be  useful  for  many  applications.  In  this  memo, 
a  closed-form  solution  for  the  phase  is  derived  as  a  function  of  the  carrier-to-noise  ratio 
(CNR).  For  the  envelope  pdf,  closed-form  solutions  are  obtained  for  different  ranges  of 
CNR  values,  which  correspond  to  the  cases  of  low  and  high  CNR,  respectively. 
Computer  simulations  are  carried  out  to  verify  the  results. 


2.  SIGNAL  AND  NOISE  MODELS 

The  signal  plus  noise,  e(t),  is 


(1)  e(t)  =  s(t)  +n(t). 

The  signal  s(t)  is  modeled  as 

(2)  s(t)  —  A[\  +mg(t)\cos2nfcft 

where  A  is  the  carrier  amplitude,  m,  0  <  m  <  1,  denotes  the  modulation  index,  g(t)  the 
modulating  function,  fc  the  carrier  frequency,  and  t  the  time  variable.  The  modulating 
function  is  assumed  to  be  zero-mean  Gaussian,  with  a  variance  of  a2n.  The  samples  of 
speech  signals  can  be  modeled  by  non-stationary  Gaussian  processes  [7]  while  the 
random  amplitude  in  radar  and  sonar  applications  is  usually  described  by  stationary  or 
nonstationary  zero-mean  Gaussian  processes  [6].  The  narrow-band  noise,  n(t),  is 
modeled  as  [1] 

(3)  n(t )  =  x(t)  cos  2nfct  —y(t)  sm2nfct, 

where  the  quadrature  components  x{t)  and  y(l)  are  independent,  identically  distributed 
zero-mean  Gaussian  processes  of  variance  a2.  They  arc  statistically  independent  of 

s(t)- 
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3.  PROBABILITY  DENSITY  DERIVATIONS 


The  signal  plus  noise  can  be  written  as 

(4)  e(t)  =  r(t)cos[2nfct+§(t)], 
where  the  envelope,  r(t),  and  the  phase,  <| )(?),  arc  given  by 

(5)  r(t)  =  \JZJ{t)+y2{t)  and  <>(?)  =  arctan||y. 

In  (5),  E,(t)  =  A  +Amg(t)  +x(t )  =  A  +  z(t).  From  the  assumptions  regarding  the 
Gaussian  nature  and  the  mutual  independence  of  g(t),  x(t)  and  y(t),  z{t )  is  zero-mean 
Gaussian  with  variance 

(6)  a:  =  A2m2o2n  +  <r  =  (2pm2a2  +  l)a2 

and  is  independent  of  y(t).  In  (6),  p  =  A2/(2a2)  is  the  CNR.  For  ease  of  exposition,  c 
and  a  are  introduced  to  be 

(7)  c  =  2pm2 a2  and  a  =  c  +  1, 

and  the  time  variable  is  henceforth  omitted.  Since  z  and  y  arc  independent,  the  joint  pdf 
of  ^  and  y,  /;(c,,  y),  can  be  written  as 

(t;-A)2  y2  - 
2  ao2  2  a2 

The  joint  pdf  of  the  envelope  and  the  phase,  p(r,ty),  is  obtained  from  the  relation 
p(r,§)drd§  =  p(t„y)d^dy.  Using  the  fact  that  q  =  rcosp.  y  =  rsin  (j  and 
dqdy  —  rdrdty,  one  obtains 


(8) 


p{t»y)  = 


2zz^J~d(5‘ 


■exp 


(9) 


p(rA) 


(1  +csin2(|))r2  —  2rAcos<|) +A2 
2  ao2 


3.1  Phase  probability  density 

The  marginal  pdf  of  the  phase,  p(§),  is  obtained  by  integrating  p(r.  (j>)  over  r.  Upon 
completing  square  in  r  in  the  numerator  of  the  exponent  in  (9),  we  have 


(10)  p{§)  = 


2%^/ac- 


■  exp 


A2nsin2(|) 

poo 

/  rpvri 

i 

<N 

1 

2fld2(l  +csin2(|))_ 

1  /  CAU 

Jo 

- 1 

<N 

0 

<N 

_ i 

dr. 


where 


2 
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(11) 


,  Acostb  .2  ac2 

b  = - j—  and  o  = - . 

1  +  c  sin-  (|)  1  +  c  sin  <|) 


The  can  be  further  written  as 


(12) 


p(4>)  = 


y/cie  + 


Via1 


1  +  erf  I 


(vl?)]  “p  {w) } 


271  (1  +csin2(j)) 


where  erf()  denotes  the  error  function  defined  by  [8] 


The  phase  pdf  has  therefore  been  expressed  in  terms  of  the  CNR  and  the  AM  signal 
parameters  m  and  a2 .  It  is  an  even  function  of  (|).  A  plot  of  p(( j>)  vs.  (f)  is  shown  in 
Figure  1  for  five  values  of  the  CNR  and  for  m  =  0.5,  a2  =  0.04.  This  value  of  a2  is 
sufficiently  small  that  most  of  the  modulating  function  fluctuations  satisfy  the  usual 
restriction  that  |g(?)|  <  1.  If  m  =  0,  then  c  =  0,  a  =  1,  and  (12)  and  (14)  reduce  to  the 
phase  pdf  for  a  sinusoidal  signal,  as  given  in,  for  example,  [1]  [p.  140,  equation 
(4-219)]. 

3.2  Envelope  probability  density 

The  marginal  pdf  of  the  envelope,  p(r),  is  obtained  by  integrating  (9)  over  (]): 
r  r2n 

(15)  PW=2 “p 

An  exact  evaluation  of  this  integral  appeal's  to  be  difficult.  Instead,  we  obtain 
approximate  expressions  for  p(r)  that  are  valid,  for  r  -C  i'o  and  for  r  >•  ip,  where  ip  is 
defined  below.  To  proceed,  we  factor  out  the  large-p  behavior  of  p{r).  When  p  I , 
rxiA  +  z.  It  follows  that  p{r)  can  be  approximated  by  an  Gaussian  function 


(1  +csin-(|))r-  —  2rAcos(|) +A“ 
2  ao2 


d§. 


(16) 


exp[ 


(r~A)2 
2 ac2  J’ 
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where  the  constant  terms  arc  ignored.  Then,  we  write  (15)  as 


(17)  p{r) 


( r  —  A )2  rn  f  cr[rsin2(|) +  2ro(l  —  cost]))] 
2 ao2  _  Jo  CXp|  2 ao2 


d( |), 


where 

(18) 


A  _  o 
c  V2p m2a2  ’ 


and  the  integral  in  (15)  over  the  interval  (0,2ji)  has  been  written  as  twice  the  integral 
over  (0, 7t).  For  r  -C  n>,  the  contribution  of  r sin2  4*  in  l*ie  integrand  in  (17)  can  be 
neglected  compared  with  that  of  2r<) ( I  —  cost]))  for  (|)  in  (0,7t).  The  envelope  pdf  then 
becomes 


(19) 


P(r ) 


r2  +  A2  \ 
2 ac2  ) 


k 


7 


where  7o()  denotes  the  modified  Bessel  function  of  the  first  kind  of  order  zero  [8].  It  is 
noted  that  if  r  and  A  are  scaled  by  i/a,  then  (19)  is  the  Rice  distribution  [1]  for  the 
scaled  variables. 

For  p  <C  1  and/or  small  signal  parameters,  c<  1,  and  the  term  c  sin 2  c)>  in  the  integrand 
in  (15)  can  be  neglected  compared  with  1.  The  resulting  p{r)  is  identical  to  that  in  (19). 

In  terms  of  the  normalized  envelope,  r'  =  r/a,  p,  and  the  signal  parameters,  (19) 
becomes 


(20) 


p(r)  ss 


r  <C  >'{)■ 


For  r  Co,  the  term  2ro(  1  —  cost]))  in  the  integrand  in  (17)  can  be  neglected  compared 
with  rsin2(|)  only  if  (|)  is  in  the  interval  (0, tc/2).  For  (f)  in  the  interval  (71/2,71),  as  it 
approaches  it,  rsin2 <])  eventually  becomes  smaller  than  2r<)(  1  —  cost])),  even  if  r  r o- 

This  shows  the  need  to  split  the  integral  into  one  over  (0,Jt/2)  and  another  over 
(7t/2,7t)  and  approximate  each  differently.  The  integral  over  (O.Jt/2),  after  neglecting 
2ro(  1  —  cost])),  is  evaluated  to  be 


The  integral  over  (jt/2,Jt)  becomes,  after  a  simple  change  of  integration  variable  and  an 
elementary  manipulation. 
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(22) 


(  2 Ar\  [n/2  f  cr[rsin2(|)  —  2ro(l  —  cos(|))]  ) 

exp(-^Ji,  exp{ - w - )«■ 

In  this  form,  the  term  2ro(  1  —  cost]))  can  now  be  neglected  compared  with  rsin2 p,  and 
the  resulting  integral  is  the  same  as  that  in  (21).  The  exponential  multiplying  the 
integral  in  (22)  is  small  compared  with  1.  To  see  this,  it  is  noted  that,  because  of  the 
Gaussian  in  (17),  p(r)  is  appreciable  only  for  r  in  the  neighborhood  of  A.  The  exponent 
2 Arj ( aa 2)  is  therefore  ~  2 A2 /(aa2)  =  4p /a.  It  is  shown  in  the  next  paragraph  that 
r  ?'o  corresponds  to  large  p,  for  which  a  ~  c.  Hence  the  exponent  is  2/{m2a;n).  In 
order  for  the  modulating  function  to  satisfy  |g(?)|  <  1  well,  the  upper  bound  on  a2 
should  be  about  0.1.  Thus  for  m  varying  between  0.1  and  1,  the  exponent  varies  from 
about  2000  to  about  20,  showing  that  the  result  of  (22)  can  be  neglected  compared  with 
that  in  (21).  Substituting  the  latter  result  in  (17),  one  obtains 


(23) 


2  aa2 


For  r  ~  A,  the  argument  of  Iq  in  (23)  is  ~  p/2  1.  Replacing  Iq  by  its  asymptotic 

behavior,  7o(x)  ~  ex /y/2nx  for  x^$>  l  [8],  yields  a  p(r)  that  is  Gaussian,  in  agreement 
with  one’s  expectation.  In  terms  of  r1,  p,  and  the  signal  parameters,  p(r)  can  be  written 
as 


(24) 


p(r)  « 


(r'-VZp)2' 

2  a 


exp 


r  >  r0. 


The  conditions  r  -C  i'o  and  r  2>  r o  can  be  considered  equivalent  to  low  and  high  CNR, 
respectively.  To  see  this,  the  maximum  envelope,  rmax,  beyond  which  p(r )  is  negligible 
is  first  estimated.  For  small  CNR,  p(r )  is  approximately  Rayleigh,  with  the  peak 
occurring  at  r  «  a.  One  may  choose  rmax  to  be  a  multiple  of  a,  say,  5a.  For  large  CNR, 
p(r)  is  approximately  Gaussian,  with  the  peak  occurring  at  r  ~  A  and  a  standard 
deviation  of  yCa.  One  may  estimate  rmax  to  be  the  sum  of  A  and  a  multiple,  again  say 
five,  of  the  standard  deviation,  resulting  in 


(25)  rmax  =A  +  5y/ac=  (\/2p  +  5V«)^- 

Since  the  high-CNR  estimate  exceeds  the  low-CNR  estimate,  (25)  is  taken  to  be  the 
maximum  envelope  for  all  CNRs.  A  comparison  of  rmax  with  ro  is  shown  in  Table  1  for 
four  values  of  the  CNR,  m  =  0.5,  and  a2  =  0.04.  In  Figure  2,  p(r)  calculated  from  the 
approximate  expressions  arc  compared  with  that  obtained  by  numerical  integration,  for 
the  same  four  CNRs  and  signal  parameters.  In  the  absence  of  any  expression  to  connect 
the  r  <  ro  and  r  3>  ro  results,  p(r)  is  calculated  from  (20)  for  r  <  ro  and  from  (24)  for 
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Table  1:  Comparison  of  f-max  with  rg ;  m  =  0.5,  o~  =  0.04. 


CNR  (dB) 

7'max  /  Ct 

r0/a 

0 

6.46 

70.71 

10 

9.95 

22.36 

20 

22.80 

7.07 

30 

67.63 

2.24 

r  >  ?•().  This  results,  in  general,  in  a  discontinuity  at  r  =  r o-  The  dashed  curves  in  Figure 
2  show  the  pdfs  calculated  in  this  way.  The  solid  curves  correspond  to  evaluating  the 
integral  in  (17)  numerically.  For  CNR  of  0  dB,  the  entire  envelope  range  is  much 
smaller  than  ro,  and  (20)  agrees  well  with  the  numerically-integrated  result.  For  CNR 
of  30  dB,  the  envelope  range  in  which  p(r )  is  significant  is  much  beyond  ro,  and  (24) 
gives  a  good  description  of  the  pdf.  For  CNR  of  10  dB,  although  rmax  <  ro,  it  is  a  much 
larger  fraction  of  r o  than  is  the  case  for  CNR  of  0  dB.  The  agreement  is  less  satisfactory, 
as  one  would  expect.  For  CNR  of  20  dB,  ro  <  rmax,  but  rmax/r o  is  much  smaller  than  is 
the  case  for  CNR  of  30  dB.  The  agreement  is  again  less  than  satisfactory.  An  analysis 
of  the  fractional  error  between  the  approximate  pdf  and  the  numerically-integrated  one 
confirms  that  the  error  increases  as  the  envelope  value  approaches  r o  from  either  side. 


4.  CONCLUSIONS 


In  this  memo,  we  have  derived  a  closed-form  pdf  for  the  phase  of  amplitude  modulated 
signals  as  a  function  of  the  CNR  and  the  signal  parameters.  For  the  envelope  pdf, 
approximate  expressions  are  obtained  for  small  and  large  CNRs.  Using  computer 
simulations,  we  showed  that  the  numerical  results  agreed  well  with  the  theoretical 
derivations.  Future  efforts  will  be  focused  on  applying  the  results  for  identification  of 
radar  and  communication  signals  for  EW  purposes. 
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